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Doping change and distortion effect on the double-exchange ferromagnetism are studied within a 
simplified double- exchange model. The presence of distortion is modelled by introducing the Falicov- 
Kimball interaction between itinerant electrons and classical variables. By employing the dynamical 
mean-field theory the charge and spin susceptibility are exactly calculated, ft is found that there is 
a competition between the double- exchange induced ferromagnetism and disorder-order transition. 
At low temperature various long-range order phases such as charge ordered and segregated phases 
coexist with ferromagnetism depending on doping and distortion. A rich phase diagram is obtained. 
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I. INTRODUCTION 

The discovery of colossal magnetoresistance in doped 
manganitesi^ has renewed interest in the ferromagnetism 
induced by the double-exchange (DE) mechanism^ The 
main feature of the DE mechanism is a cooperative ef- 
fect where electron hoping favors ferromagnetic (FM) 
ordering of localized spins via the FM Hund coupling, 
and vice versa, the presence of the FM order facilitates 
the electron hoping. The occurrence of metallic FM 
state in doped manganites Ri_^A^Mn03 (where R is a 
trivalent rare-earth element and A is a divalent alkaline 
ion) was qualitatively explained by the DE mechanism. 4 
The physically relevant electrons in these compounds are 
those from Mn 3d levels which are split by the cubic crys- 
talline field into triply degenerate t 2g levels and doubly 
degenerate e g levels. Electrons of e g levels are able to 
hop between Mn sites and form the conduction band, 
while electrons of t2 g levels are localized. Conduction 
electrons and localized spins are correlated by the DE 
mechanism which leads to the appearance of the metal- 
lic FM phase. The DE model became the starting point 
toward comprehensive understanding of the properties of 
doped manganites. 

Experiments have observed in doped manganites a very 
rich phase diagram, which involves phases with spin, 
charge and orbital ordersi For the undoped case (x = 0) 
all Mn ions are Mn 3+ and are expected to induce a Jahn- 
Teller (JT) distortion. For the other extreme doping 
(x = 1) all Mn ions are Mn 4+ and do not couple to 
the JT distortion. In the regime of intermediate dop- 
ing, two valence ions Mn 3+ and Mn 4+ are simultane- 
ously present. The presence of two valence ions may 
lead to a static mixed valence Mn 3+ /Mn 4+ configura- 
tion, in particular, to an alternation charge-ordered (CO) 
state of Mn 3+ /Mn 4+ ions for appropriate dopings. This 
is the conventional view, for which there are abundant 
experimental and theoretical supportsi^ In particular, 
recently a CO-FM state has been observed^ However, 
there are several experiments which challenge the con- 
ventional view. Several x-ray absorptionSi&fi and neu- 
tron diffraction studies^ revealed pictures that do not 
match with the static mixture of Mn 3+ and Mn 4+ ions. 



One suggests that all Mn ions have the same valence 
and result into the Zener-polaron stateiii*^ However, 
very recent experiment observed the presence of two 
types of Mn sites with different local geometric struc- 
tures. One of the types of Mn sites is surrounded by 
a tetragonal-distorted oxygen octahedron, whereas the 
other has a regular octahedral environment. As a re- 
sult a charge segregation state was deduced. With the 
motivation of the experimental observations^ we model 
the presence of the two types of Mn sites by incorporat- 
ing the Falicov-Kimball (FK) modeli 3 . The FK model 
was initially introduced as a statistical model for metal- 
insulator transition^ 3 " Later it was also applied to va- 
lence change transitions in intermetallic compounds^ 4 ^ 
Within the FK model the presence of two types of Mn 
sites is mapped to a classical variable which only ac- 
cept two values (for instance, 1 and 0). The energy 
difference of these sites is mapped into the interaction 
strength of the FK model. Indeed, the sites surrounded 
by tetragonal-distorted octahedron have induced the JT 
distortion. As a consequence the energy levels of the 
distorted Mn sites are split. The FK model can de- 
scribe a charge ordered phase as well as a charge seg- 
regated phased In particular, the model can exhibit the 
checkerboard CO state in appropriate conditions. The 
checkerboard CO state is truly a mixed-valence state. 
The segregated state is a phase-separated mixture of two 
full uniform configurations>iZii2ii2i2£ In such the way, at 
low temperature the FK model could establish various 
phases with different charge configurations which may 
correspond to the experimental observations. 8 

However, the FK model alone cannot describe the DE 
induced FM state which was also observed in doped man- 
ganites. Therefore we incorporate the FK model into 
the DE model in order to study both the charge ordered 
phases and the ferromagnetism upon doping. The com- 
bined model has previously been considered in the con- 
text of order-disorder change of the A-site substitution. 21 
In the previous studjiSi only the checkerboard CO and 
FM state are considered in the limit of infinite value 
of the Hund coupling. In this paper we study all pos- 
sible ordered phases of the combined model in whole 
range of doping and interaction. In order to detect the 
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phase transition we study the static charge and spin re- 
sponse of system within the dynamical-mean field theory 
(DMFT)iSa The DMFT has been widely used for investi- 
gating strongly correlated electron systems. Within the 
DMFT the static charge and spin correlation function are 
calculated explicitly. We find that the system exhibits a 
rich phase diagram which includes various charge ordered 
phases coexisting with ferromagnetism. In particular, the 
checkerboard CO state or the segregated state can coexist 
with the FM state. The combined model can also serve 
as a model for studying the problem of order-disorder A- 
site substitutional or the problem of orbital ordering in 
doped manganitesi 23 i 24 

The present paper is organized as follows: In Sec. II we 
present the combined FK and DE model and its DMFT 
solutions. In Sec. Ill the charge and spin correlation 
functions are calculated explicitly. The numerical results 
and discussions are presented in Sec. IV. The final section 
is conclusion and remark. 



II. MODEL AND DYNAMICAL MEAN-FIELD 
THEORY 

The combined FK and DE model in our study is de- 
scribed by the following Hamiltonian 

<r/>,<7 ia i 

+U'^n ic ,Wi + Ew^^Wi, (1) 



where c\ a (c ia ) is the creation (annihilation) operator for 
an itinerant electron with spin a at lattice site i. The 
first term in Hamiltonian Q represents the hoping of 
itinerant electrons between the nearest neighbor sites, t 
is the hoping integral and is scaled with the spatial di- 
mension d as t = i*/(2v / d)^ In the following we will 
take t* = 1 as the unit of energy. Sf is the z com- 
ponent of localized spin at lattice site i. For simplic- 
ity, it takes two values —1,1. s\ = (ct-Cj-f — cj|C^)/2 
is the z component of itinerant electron spin. Wi is a 
classical variable that assumes the value 1(0) if site i is 
surrounded by distorted (regular) octahedron. U is the 
interaction strength and is mapped into the difference 
in the level energy of these sites. The expectation value 
p w = J2i( w i) /N, {N is the number of lattice sites), corre- 
sponds to the concentration of distorted sites. The chem- 
ical potential \x controls the doping n e = ^2 ia .{njq-} / f N , 
while E w controls the fraction of distorted sites. The con- 
dition n e + p w — 1 is used to determine E w for each dop- 
ing n e . The first three terms in Hamiltonian describe 
a simplified DE model which contains only the Ising-type 
interaction between the itinerant and localized electron 
spins. The simplification does not allow any spin-flip 
processes, which can be important at low temperature 
where spin-wave excitations may govern the thermody- 
namics of the system. However, in the DE processes spins 



of itinerant electrons align ferromagnetically with the lo- 
calized spins, hence the Ising part of the Hund coupling 
plays a dominant role. Within the DMFT the simplified 
DE model is equivalent to the DE model with classical 
localized spins in the disordered paramagnetic phase, 4 
The transport quantities calculated within the simplified 
model capture essential features of the full DE model. 27 
The simplified DE model has also previously been used in 
the study of doped manganitesi 2 ^ The last two terms in 
Hamiltonian take into account the energy difference 
of Mn sites. They together with the hoping term form 
the FK models Several authors have also constructed 
the combined model to study the properties of mangan- 
ites in different contexts and regimesi 21 * 2 ^ 2 ^* 2 ^ Ferrari 
et al. used the combined model to study the metallic 
FM phase of the two orbital DE models Recently, Ra- 
makrishnan et al. basically used the combined model to 
construct a two band model of localized polaronic and 
broad band states^ They used the DMFT to calculate 
transport quantities and explained the metal insulator 
transition and the colossal magnetoresistance in doped 
manganites^ 

We solve the combined model JT} by the DMFTi 2 ^ The 
DMFT is based on the infinite dimension limit. In the 
infinite dimension limit the self-energy is pure local and 
does not depend on momentum. The Green function of 
itinerant electrons with spin a satisfies the Dyson equa- 
tion 
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e(k) + fi- Y, a (iuj n y 



(2) 



where e(k) = — 2t J2i=i d 00 ^*) * s * ne dispersion of free 
itinerant electrons on a hypercubic lattice, and E^iu^) 
is the self energy which depends only on frequency. The 
self energy is determined by solving an effective single- 
site problem. The action for this effective problem is 

S c g=- dr rf'r / y^c^(r)g~ 1 (r - r')c CT (r') 
Jo Jo a 

- / dTJ2[J H S z a + /i - Uw]cI(t)c*(t) + f3E w w, (3) 

where Q a (j — t') is the Green function of the effective 
medium. It plays as the bare Green function of the effec- 
tive problem. The local Green function also satisfies the 
Dyson equation 



(4) 



where Q a (iuj n ) is the Fourier transform of Q a (r). The 
local Green function G <T (iw n ) of the effective single-site 
problem is solely determined by the partition function 



SQ~ (iuj n ) ' 



(5) 



where Z e g is the partition function of the effective prob- 
lem ©. The self-consistent condition of the DMFT re- 
quires that the local Green function G a (ioj n ) obtained 
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within the effective problem must coincide with the local 
Green function of the original lattice, i.e., 



k 

dep(e) 



1 



iu) n - e + fi- S ff (iw„) : 



(6) 



where p(e) is the density of state (DOS) of noninter- 
acting itinerant electrons. In the infinite dimension 
limit of hypercubic lattices it has the form p(e) — 
exp(-e 2 /(i*) 2 )/\/7rf*- Eqs. Q, © and © form the 
self-consistent equations for determining the self-energy, 
and hence, also the Green function of the original lat- 
tice. Within the effective single-site problem, the parti- 
tion function is 



Z, 



off 



Tr / DclDc^e 



(7) 



where the trace is taken over S z and w. This partition 
function can be calculated exactly. It is similar to solve 
the FK model within the DMFT.22 We obtain 

Z cS = 2 ex P [ ~ P E ^a + 

Zuiiuin) + asJn — all 



a = 0,l s = ±l 



5> 



,(8) 



where Z a (iu> n ) = Q a l (iuj n ). Using Eq. |J5J| we obtain the 
local Green function 



G a (iuj n ) 



^ Z a {iuj n ) + ijsJh — all ' 



(9) 



where the weight factors W as are 



loff 



■ exp 



-j3E w a- 



Z a {iuj n ) + ijsJh — otU 



(10) 

with a — 0, 1 and s = ±1. Note that the weight factors 
W as are not simply a number. They are functionals of 
the local Green function. This is an important feature 
of the DMFT that gives nontrivial contributions to the 
response functions of the system. 29 In the paramagnetic 
phase Z^(iui n ) — Zi(iui n ), hence W as = W a ,- S that leads 
the Green function (J5J) and the self energy are indepen- 
dent of spin indeces, as expected. The value of E w is 
adjusted that the concentration p w fulfills n e + p w = 1 
for each doping n e . One can show that 



fh, 



s=±l 



(11) 



We use this equation to adjust the value of E w . So far, 
we have obtained closed system of equations for deter- 
mining the Green function of the system. The system of 
equations can be solved numerically by iterations:- 6 ' 



III. INSTABILITY OF HOMOGENEOUS 
PARAMAGNETIC PHASE 

In order to detect the charge and spin ordered states 
which are established at low temperature we study the 
static charge and spin correlation function of itinerant 
electrons in disordered paramagnetic phase. The signal 
of a phase transition is a divergence of these correlation 
functions at a certain momentum. The charge (c) and 
spin (s) correlation function are defined as 

X c(s) (hj) = ((Sm^SnnXSnj^Snji)), (12) 

where brn a — — (tw)- These correlation functions 
can be expressed as 



(13) 



where XoV {h j) = (Srii a dnj a i), and £ CT = 1 for the charge 
correlation function and £ CT = a for the spin correlation 
function. In order to calculate the static correlation func- 
tions, one has to introduce external fields h.^ which cou- 
ple to riic into the Hamiltonian. The correlation func- 
tions Xerer'ihj) are obtained by differentiating the Green 
function with respect to the external fields and then tak- 
ing the zero limit of the fields^ i.e., 



Xaa'(i,j) = ~T 2 Y 



dG iija (iuj n ) 



dh 



16.29 



(14) 



{h}=0 



Following the standard technique ^imtuk one can express 
the charge and spin correlation function in the terms of 
the charge and spin susceptibility x^'fa) 



X c(8) (q) = -T 2 ^x c(s) (q,^n), 



(15) 



where x ^ (q) is the static charge (spin) correlation func- 
tion in momentum space. From the definition of charge 
and spin correlation function (|13|) and relation 114(1 . we 
obtain 

X c(s) (q,«w„) = 2xo(q,«w„) + 

Xo(q,^„)i ]T |^Nr^X c(s) (q^),(16) 

vaa' 

where Xo(q, iu n ) = G ff (k + q, MJ„)G CT (k, iu n ) is the 
bare particle-hole susceptibility. Here we have used the 
fact that in the paramagnetic phase 

dEi(iu) n )/dG^(ioj n ) = d^i(iu n )/dGi(iu} n ), (17) 



dY,j(iu} n )/dGi(iuj n ) = dEi(iu n ) / 'dG^(ioj n ). (18) 

From Eqs. Q and ©, we obtain the self-energy T, a (ioj n ) 
as a functional of G a (ioj n ) and W as , hence its full deriva- 
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tives in Eq. (|16f) are expressed through its partial deriva- 
tives. We obtain 

iG a i{iuj v ) nv \dG a >(iuj n )J 



dG a i(iu> u ) 

'dT, a (icj n 



E 



(iu) n )Jw 
SW as 



dW as J G,w-^5G a '(iuj n 



(19) 



where Was means all weight factors W except W as . Sub- 
stituting l|19|) into Eq. we arrive at 



X c(s) (q,^«) 



2 1 \p / dT, a {iuj n 

2 2^V dW as Jg.w- 



^67as(q) 



_i 1 f^<y( iuj n) 



[ Xo (q,iuJ n )} 1 - - ( 



dG a ) w 

(20) 



where the matrix elements of 7(q) are 



The functional derivative of W as in the above equation 
can explicitly be expressed through the derivatives of 
W as with respect to Z a (iuj n ) and the partial derivatives 
of the self energy Y, a {iuj n )w& Substituting (|2T)jl into J2U 
and performing some algebraic calculations we obtain the 
following matrix equation 



Q(q)7(q) = P(i), 



(22) 



where the matrixes Q(q), -P(q) have the following ele- 
ments 



J 



Qas,a' s' (q) — $as : a' s' H~ 

i? QSi(T (^„)(G 2 (zc n )(|§^4) -l) 
P„(q) = ^W^U ■ ( 24 ) 



Here we have introduced the following notations 



j? C \ dWas c c 

aZ a \w n ) 

OW as 

and ?? n (q) = — G~ l (itJ n ) + G(icj n )x^ 1 (q, iw n ). In Eqs. (|2*5|) - (|2~lf> the spin indices of the Green function and the self 
energy are omitted since they are in the paramagnetic phase. The derivatives in Eqs. (|23[) - l|24|) can be calculated 
explicitly. Straightforward calculations give 



dTAwn) _ 1 

dW as (Z(iui n ) + asJ H - aU)A s (iuj n ) ' 

dT,(iuj n ) _ A G (iuj n ) 
dG(iuj n ) Aziiujn)' 
dW as W as 



(25) 
(26) 



where 



dZ a (iuj n ) Z{iui n ) + usJh — ctU 

W^ 

(Z(iu> n ) + ijsJh — cell) 



W as G(iu n ), (27) 



2 ' 



A G (uu n ) = -VV^ t-t- G [ Z ( lLUn ) + (7S J H -aU -G-\nu n )]. 

' (Z(iLu n ) + (TsJ H - aljy 



OS 



In such the way, Eqs. (19)-(24) fully determine the susceptibilities, once the self-consistent equations of the DMFT are 
solved. The charge or spin correlation function diverges whenever 7~(q) diverges, which happens when the determinant 
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of Q(q) vanishes. The divergence indicates the instability of the disordered paramagnetic state. The q dependence 
of the susceptibilities comes entirely from r? n (q), and hence from the bare susceptibility x (q, iw„). Within the 
DMFTiS^ in the infinite dimension limit all of the q dependence of the bare susceptibility can be summarized in a 
single parameter X = J2a=i cos Qa/d. One can show thaUS^ 

Xo(q,tuJ») = Xo{X,iU3 n ) = — / de- — — , 

where .F^z) = J dep(e)/(z — e) is used to denote the Hilbert transform. Now the instability of the disordered 
paramagnetic phase happens whenever the determinant of Q(X) vanishes at a certain value X. In particular, X = — 1 
corresponds to the checkerboard zone-boundary point q = (tt, tt, ...,7t) and the corresponding instability leads either 
to the charge checkerboard phase or to the antiferromagnetic phase at low temperature. X = 1 corresponds to the 
uniform zone center point q = and the corresponding instability leads to a charge segregation phase or FM phase 
at low temperature. 



IV. NUMERICAL RESULTS 

First, we consider the magnetic instability. In this case 
we calculate the spin correlation function as a function 
of X and temperature T. A divergence of the spin cor- 
relation function indicates a magnetic instability. For 
most values of Jh, U and n e the spin correlation func- 
tion diverges only at X = 1. In Fig. 1 we plot the typical 
behavior of the spin correlation function. The divergence 
of the spin correlation function at X = 1 indicates the 
FM stability. This means the FM state is established at 
low temperature. However for small values of Jh and n e 
closed to 1, the spin correlation function also diverges at 
X = — lpiiliiL This divergence indicates the stability of 
the antiferromagnetic phase at low temperature. In this 
paper we only consider the FM phase induced by the DE 
mechanism and its coexistence with CO phases. Thus in 
the rest of paper we consider the FM stability only. In 
Fig. 2 we present the FM transition temperature Tp as 
a function of U and n e for various values of Jh- Tp is 
determined from the vanishing condition of the determi- 
nant of Q(X) at X — 1. Fig. 2(a) shows that the critical 



temperature Tp decreases as increasing U and increases 
as increasing Jh- One expects in the limit Jh — > 00 the 
FM transition temperature reaches its maximum value 
for fixed U. For fixed Jh the FM transition temperature 
is maximum if there is no JT distortion (i.e. U = 0). 
The JT distortion splits the energy level of Mn ions, and 
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FIG. 2: The FM transition temperature Tf as a function of 
FIG. 1: The spin correlation function as a function of X at U [(a)n e = 0.5] and as a function of n e [(b)U=0.5] for various 
different temperatures (n e — 0.5, Jh = 2, U = 0.5). values of Jh- 
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FIG. 3: The charge correlation function as a function of X 
at different temperatures, (a) Jh = 0.1, (b) Jh = 0.2, (c) 
Jh = 0.4, (d) J H = 0.5. In all figures U = 1.0, n e = 0.6. 



this leads to suppress the FM transition temperature. 
This reduction of the FM transition temperature due to 
distortion is gradually significant already at intermedi- 
ate values of U . In Figure 2(b) we also present the FM 
transition temperature as a function of doping n e . In the 
limit Jh — > oo, the FM transition temperature is max- 
imum at half doping n e = 0.5^i However, for finite Jh 
its maximum shifts away from the half doping, to lower 
doping region. 

Next, we consider the charge ordering instability. In 
this case, we study the divergence of the charge corre- 
lation function in the homogeneous paramagnetic phase. 
In Fig. 3 we plot the typical behaviors of the charge cor- 
relation function. They show that the charge correlation 
function may diverge at X = 1, X = —1 or at an inter- 
mediate value — 1 < X < 1. The divergence at X = — 1 
indicates the checkerboard charge ordered state estab- 
lished at low temperature, while the divergence at X = 1 
indicates the segregated state established at low temper- 
ature. The divergence at an intermediate value of X in- 
dicates the charge ordered phase being incommensurate 
at low temperature. The charge ordering critical tem- 
perature T c is determined from the vanishing condition 
of the determinant of Q(X). However, one notices that 
for fixed values of Jh, U and n e the determinant of Q(X) 
may vanish at different X and temperature T. Hence, we 
obtain the critical temperature T C (X) as a function of 
X. However, this does not indicate the stability of many 
charge ordered phases at low temperature. For certain 
values of Jh, U and n e there is only one charge ordering 
transition which happens at the maximum temperature 
T c among T C (X). Below this temperature T c , although 
the charge correlation function may still diverge at other 
values of X, the divergence does not indicate a charge 
ordering, because the assumption of disordered phase is 
not valid anymore. Therefore the charge ordering critical 
temperature is determined not only from the vanishing 



FIG. 4: The dependence of T c on X in case of U = 1.0, 
n e = 0.6 with various values of Jh- 



condition of the determinant of Q(X), but also from the 
maximum condition^ 

{ x} cV 1 

The corresponding value of X at which T C (X) is maxi- 
mum determines the charge arrangement of the low tem- 
perature phase. In Fig. 4 we plot function T C (X) for 
different values of Jh- It shows that for fixed Jh, U, n e 
we always find a unique maximum T c . For small values 
of Jh wc obtain the segregated state at low temperature. 
For large values of Jh the checkerboard ordered phase is 
observed. In an intermediate regime we also find an in- 
commensurate charge ordering phase transition. In Fig. 
5 we plot the charge ordering critical temperature as a 
function of U for various values of Jh- It shows that 
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FIG. 5: The critical temperature of charge ordered phase 
transition as a function of U (n e = 0.5). The filled (opened) 
symbols and the dotted lines correspond to the checker- 
board (incommensurate) charge ordered phase and segregated 
phase, respectively. 
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FIG. 6: The charge ordering critical temperature T c as a func- 
tion of Jh (U — 1.0). The filled (opened) symbols and the 
dotted lines correspond to the checkerboard (incommensu- 
rate) charge ordered phase and the segregated phase, respec- 
tively. 



the checkerboard ordered phase is established for large 
values of Jh- As increasing U, the critical temperature 
first increases, reaches a maximum value, and then de- 
creases. The behavior of the critical temperature is sim- 
ilar to the one obtained in the FK model^ For small 
values of Jh the checkerboard CO phase disappears and 
the segregated phase or incommensurate CO phase may 
be established depending on the value of U, as shown in 
Fig. 5. Figure. 5 also shows that the segregated phase 
is established at large values of U, while the incommen- 
surate phase is established at smaller values of U. One 
notices that the segregated phase detected from the di- 
vergence of the charge correlation function constitutes a 
continuous phase transition. However, the phase transi- 
tion is indeed first orderii&S It can be shown by consid- 
ering the free energy and using a Maxwell construction 
at low temperatureiiSiS Usually, the critical temperature 
of the first order phase transition is higher than the one 
obtained from the divergence of the charge correlation 
functioni 16 ! 32 

In Fig. 6 we plot the critical temperature T c as a 
function of Jh for various doping n e . It shows that the 
checkerboard charge ordered phase is established at large 
values of Jh and disappears at small Jh ■ Its critical tem- 
perature increases as increasing Jh- In contrast, the seg- 
regated phase is established at small values of Jh and its 
critical temperature decreases as increasing Jh- For in- 
termediate values of Jh, both the checkerboard CO and 
segregated phases disappear, and an incommensurate CO 
phase is established. In Fig. 7 we plot the critical tem- 
perature T c as a function of doping n e for various values 
of Jh ■ For large values of Jh the checkerboard CO phase 
is established, and its critical temperature reaches maxi- 
mum at half doping n e = 0.5. For intermediate values of 
Jh, the CO phase is established at large doping n e > 0.5, 
and disappears at smaller doping. Instead of the checker- 



FIG. 7: The charge ordering critical temperature T c as a func- 
tion of n e (U — 1.0). The filled (opened) symbols and the dot- 
ted lines correspond to the checkerboard (incommensurate) 
charge ordered phase and the segregated phase, respectively. 



board CO phase, an incommensurate phase appears. For 
small values of Jh only the segregated phase appears at 
small doping. At large doping the segregated phase also 
disappears. 

We summarize the above results with the phase dia- 
gram in Fig. 8 which plots the regions of stability for 
different charge ordered phases. The stability is deter- 
mined by the symmetry label X of the initial ordered 
phase as the temperature is lowered to the first instabil- 
ity at T c . Furthermore, we assume that the symmetry 
label X of the ordered phase does not change as the tem- 
perature lowered from T c to zero. Actually, the phase 
diagram is an approximation of the zero temperature 
phase diagram^ The phase boundaries may change as 
one reduces the temperature from T c to zero since the 
behaviors of incommensurate phases at low temperature 
are not able to be studied within the present approach. 
Moreover, the phase boundaries may also change if there 
are first-order phase transitions which may happen with 
the segregation phases. The phase diagram shows that 
the incommensurate phases is stabilized in buffer zones 
between the disordered and checkerboard CO phase or 
between the disordered and segregated phase. The seg- 
regated phase exists only for small values of Jh ■ For large 
values of Jh the checkerboard CO phase is stabilized. So 
far, we have obtained different charge ordered phases de- 
pending on the value of Jh, U and doping n e . On the 
other hand, the system always exhibits the FM stability. 
Although the FM transition may happen first, in the FM 
phase the charge density still remains homogeneous and 
the charge correlation function is still normal as in the 
homogeneous paramagnetic phase. Hence one can use the 
charge instability signal in the high temperature homo- 
geneous phase as the signal of charge ordering even if the 
system is already in the FM phase. Indeed, the charge 
ordering critical temperature detected from the charge 
order parameter in the CO-FM state coincides with the 
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regated phase coexisting with the FM phase is obtained 
theoretically. The phase is a phase separated mixture 
of two types of Mn sites with different local geometric 
structures. The regime of the phase separation is partic- 
ularly interesting, and possible consequences of its exis- 
tence may be relevant to the experimental observations & 
However, the segregated phase appears only for weak 
Hund coupling and strong distortion. There is also a 
possible coexistence of the two segregated phases, one of 
which is between the magnetic phases, and the other is 
between the charge ordered phases. However, we leave 
the problem for further study. 
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FIG. 8: Charge ordering phase diagrams: (a) (7 = 1; (b) 
Jh = 1- The shorthand CO denotes the checkerboard charge 
ordered phase; IC, the incommensurate phase; SP, the segre- 
gated phase, and CD, the charge disordered phase. 



one calculated from the instability of the charge correla- 
tion function in the homogeneous paramagnetic phase, — 
In such the way, the FM state coexists with different 
charge ordered phases for appropriate doping and distor- 
tion. One notices that in the previous studies^i*2i2iS a 
rich phase diagram which includes spin, charge and or- 
bital ordered phases was also obtained. There is also a 
coexistence of the FM phase with checkerboard CO phase 
due to the Jahn- Teller phononsi 31 i 33 i 34 i 35 In particular, 
the ferromagnetic CO phase is stabilized for large Jahn- 
Teller coupling and infinite Hund coupling. The con- 
ditions are in an agreement with the phase diagram in 
Fig. 8. However, in the previous studies^ 3 ^* 3 ^ 3 ^ only a 
phase separation between different magnetic phases was 
considered. A separation of charge ordered phases has 
not been addressed. In the present paper, a charge seg- 



In the present paper we have considered the doping and 
distortion effect on the double-exchange ferromagnetism. 
By employing the DMFT we have exactly calculated the 
charge and spin correlation function. A long range or- 
der is determined from the divergence signal of the cor- 
relation functions. The obtained results show that the 
system exhibits various phases which include the FM, 
checkerboard CO, incommensurate CO and segregated 
phases. In particular, the FM phase can coexist with the 
checkerboard CO phase for large values of Jh and with 
the segregated phase for small values of Jh ■ The incom- 
mensurate phases appear in the buffer zones between the 
regions of the charge ordered phases with different sym- 
metries. It is interesting to note that experiments have 
observed both the CO-FM stated and charge segregated 
phase. 8 The phase separation between different charge or- 
dered phases is a novel regime in manganites. By includ- 
ing the distortion effect in manganites via the Falicov- 
Kimball interaction we have simulated the charge segre- 
gated phase. The phase diagrams were found to clearly 
distinguish regions with robust CO-FM correlations and 
charge phase separation. However, manganites are too 
complicated of a system to be completely described by 
this simple model. In particular, experiments have ob- 
served inhomogeneous regions with different long-range 
orders^ The study of the properties is beyond the capac- 
ity of the present method. 
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